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SUMMARY 



A theoretical investigation is undertaken to study the buckling of a 
thin^ rectangular^ homogeneous^ isotropic flat plate augmented with 
a centrally located longitudinal stiffener. The influence on overall 
plate efficiency of the use of advanced composite materials in the 
stiffener is developed and discussed. Results of the analysis are 
presented in the form of charts to demonstrate the effects of fibrous- 
composite reinforcement of the stiffener on the critical-s tress- to- 
weight ratio of the plate. 
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INTRODUCT ION 



The trend to higher performance aircraft in all speed regimes has led 
designers to seek methods of decreasing the structural weight of air- 
craft while maintaining or increasing the stress levels on the mate- 
rials involved. One of the most promising means of accomplishing this 
objective has been the introduction of advanced composite materials to 
either augment or replace conventional structural materials (for exam- 
ple^ aluminum). One type of composite material consists of the in- 
corporation of high-strength^ high-modulus filaments into a compatible 
matrix. Several material system classes currently of active interest 
include graphite epoxy^ boron epoxy^ and metal matrices having metal 
fibers. As pointed out in Reference 1^ design and structures engi- 
neers^ using composites^ have the ability to build directional prop- 
erties into the structure which^ in turn^ results in improved efficiency. 
Weight saved^ with its direct effect on the design index of strength/ 
weighty is the most attractive feature of the use of high-performance 
composites. It has not been uncommon for s tructures-materials engi- 
neers to predict weight reductions of as much as forty percent in total 
systems weight through the potential use of advanced composites in both 
primary and secondary structures. 

With advanced filamentary composites^ increased design flexibility of 
significant degree exists that is not possible with conventional metal 
structures. However^ as with many new developments^ work in the com- 
posite field has been primarily empirical^ while theory and analysis 
have lagged somewhat behind. Despite the potential of composites^ 
evidence of their use appears only in the reasonably small structural 
elements that can be designed by test (for example^ stabilizers^ wing 
tips^ slats^ and rotors). The trial and error procedure^ especially 
with materials of such high performance possibilities as composites^ 
can lead to excessive costs due to both inherent high material costs 
and eventual appearance of structural behavior neither anticipated nor 
properly estimated. One such case is presented by Durlofsky and Mayers 
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in Reference 2. Another such case exists with the use of advanced 
composites to enhance the buckling efficiency of stiffened plates 
utilized in the construction of wing^ fuselage^ and empennage struc- 
tures of aircraft. Peterson^ studied the effect of an all-composite 
stiffener in the reinforcement of a simply supported aluminum plate 
without considering the possibility of an optimized composite-metal- 
stiffener combination. His results indicate that some caution must 
be exercised in using high-performance composites in combination with 
conventional structural concepts and prediction methods. 

Earlier studies of conventional stiffened plates under compression^ 
such as those of Barbre^ Seide and Stein^^ and Gomza and Stein^ 
neglect the effects of torsional stiffness of the stiffener in es- 
tablishing buckling criteria. The justification offered is usually 
that the torsional stiffness of the stiffeners is relatively small 
when compared with the bending stiffness and that vanishingly small 
torsional stiffness of the stiffeners is consistent with the appear- 
ance of buckle nodes at stiffener locations in the chordwise direc- 
tion. Such justification^ of course^ is appropriate for the discrete 
stiffener spacing of low subsonic speed flight vehicles. For the 
higher speed ranges^ stiffeners are more closely spaced and the 
stiffener torsional stiffness becomes significant in buckling anal- 
ysis. Thus^ the introduction of composite materials into all or part 
of a stiffener can possess ramifications not evident in the use of 
conventional materials. For example^ although Young’s modulus (e) 
for a composite may be as much as four times greater than that of a 
structural metal^ the shear modulus (G) may be lower by approxi- 
mately the same factor. Therefore^, neglecting torsional rigidity 
in a composite-metal system with closely spaced stiffeners may in- 
troduce unsatisfactory errors in buckling calculations. 

It is the purpose of this investigation to (l) present the theory and 
an approximate analysis for determining the stability of a simply sup- 
ported^ flat^ rectangular orthotropic plate with a centrally located^ 
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composite-augmented^ longitudinal stiffener and (2) determine for a 
square plate the effects of varying the percentage (by volume) of 
composite material in the stiffener with stiffener torsional rigidity 
included. 
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GENERAL THEORY 



STRAIN-DISPLACEMENT RELATIONS 



The strain-displacement relations for small strain, large displacement, 
and moderate rotation behavior of an ortho tropic plate (see Figure l) 
can be obtained from Reference 7 
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With a single stiffener buckling into no more than 2 ha If -wave lengths 
and a homogeneous, isotropic base plate, however, interlaminar and 
transverse shear effects become of higher-order significance and they 
are assumed to be negligible throughout the analysis* Thus, the trans- 
verse shear angles 7 ^ and 7 ^ are set equal to zero. Further, since 
buckling is of interest, the stretching of the middle surface of the 
pi ate (and stiffener) may be neglected; thus 
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STRESS-STRAIN RELATIONS 



For an orthotropic plate in which the stresses do not exceed the elas- 
tic limit, Hooke’s law is applicable; it is expressed in the following 
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form relating normal force per unit length of cross section to strain: 
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A single longitudinal stiffener possesses the same stress-strain re- 
lationships as a one-dimensional beam; that is 
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POTENTIAL ENERGY 



The strain energy stored in an elastic plate of constant thickness due 
to bending and twisting may be written as 
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Substitution of Equations (2) and (3) into Equation (5) and subsequent 
integration over the thickness of the plate yields 
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With the use of the Dirac delta function^ the potential energy of a 
single centrally located symmetric longitudinal stiffener (see Fig- 
ure l) may be described over the same interval as 
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Substitution of Equations (2) and (k) in Equation ( 7 ) and subsequent 
integration over the stiffener height yields 
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If an axial^ inplane load^ positive in compression^ is applied to the 
edges of the plate^ the potential of this applied load (including the 
end load on the stiffener) is 
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However^ since in the absence of middle surface stretching (c^ ~ “ 

7 ^^ = 0 when z = O); it can be observed from Equation (l) that 
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Substitution of Equation (lO) in Equation ( 9 ) yields 



L b 



V = 



-Iff 



0 0 



— 2 — 2 — 

N w + N w + 2N w w 
X ,x y ,y xy ,x ,y 



N 



(A ) 

^ s ' 



V * (’ 



dxdy 



( 11 ) 



Finally^ the total potential is the sum of Equations (6)^ (8)^ ^nd (ll) 
and becomes 



L b 



U + V = 



= ;//h 



D ^ D |i + D [JL 

X __2 X y y X 



0 0 



\i p. 

X y 



w + 

.XX 

1 - 



w w 
,XX ,yy 



X' y 



1 - p p ^ 

^x^y 



2 2 
w + 2D w + 
yy xy ,xy 



(El) + (GJ) 



, XX 



xy 






- N ■ N - 2N w 
X ,x y ,y xy 



w ^ (A ) 6 

,y ^ ,x y 



dxdy 

( 12 ) 



VARIATIONAL PRINCIPLE AND GOVERNING EQUATIONS 



In accordance with the potential energy principle, the total potential, 
U + V , must have a stationary value with respect to admissible var- 
iations in the lateral deflection, w , in order for equilibrium to 
be achieved. This variation leads to the equilibrium or Euler equation 
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of the variation^ plus the applicable boundary conditions. The var- 
iation of Equation (12) leads to the following equilibrium equation: 
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The associated boundary conditions for the case of single support are 
obtained as 
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The solution of Equation (15); consistent with the boundary conditions. 
Equation (ih) leads to the deflection function, w . This, however, 
is a difficult procedure in general and is not the approach followed 
in this investigation. Rather, the direct method of the calculus of 
variations (Rayleigh-Ritz method) is used to determine the buckling 
criteria for w expressed in terms of a series of arbitrary parameters. 
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METHOD OF SOLUT ION 



With the total potential energy given by Equation (12) in terms of the 
plate deflection w ^ the direct method of the calculus of variations 
(Rayleigh-Ritz method) is used to determine the minimum buckling load 
for varying stiffener size and composite content. 

GENERAL CASE 



A deflection function satisfying the boundary conditions for a flat^ 
simply supported^ stiffened rectangular plate under axial compression 
is 
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Substitution of Equation (I 5 ) into Equation (12) (with N^ = N^^ = 0 ) 
and subsequent integration over the surface of the plate yields (see 
Appendix l) 
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The vanishing of the first variation of the total potential energy with 
respect to ^ the nondimensional arbitrary displacement parameter^ 

yields the equilibrium condition and a criterion for obtaining the 
buckling load. This procedure^ as shown in Appendix I, yields 
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HOMOGENEOUS, ISOTROPIC SQUARE PLATE 



A homogeneous^ isotropic square plate with an equivalent stiffener (to 
account for varying amounts of composite material) has been selected 
as the specific case for examination in the present investigation. A 
square plate with one centrally located longitudinal stiffener under 
these loading conditions would buckle in one of two configurations: 

(l) one half-wave in each direction (m = n = l) or, (2) two half- 
waves in each direction (m = n = 2) . Any other combinations of 

ra and n would not lead to the minimum value of the buckling co- 
efficient k^ . As noted from Equation (l?) the bending stiffness 
term (El)/Db of the stiffener is effective only during the 1,1 
mode, and conversely, the torsional stiffness term (GJ)/Db is ef- 
fective only during the 2,2 mode. 

For an approximate solution, the first term of the deflection function. 
Equation (l5); could be used. However, this would allow only for a pure 
sinusoidal half-wave deflection pattern of the plate in the chordwise 
direction. As stiffener size or stiffness is increased, rotation at the 
inflection point of the 2,2 mode is somewhat restrained, indicating the 
need for at least a two-term deflection function; that is, n = 1,3 for 
the 1,1 mode and n = 2,4 for the 2,2 mode of Equation (l?)* The 
result of carrying out this procedure (see Appendix l) is 
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for the 2^2 mode^ where D is the bending rigidity for a homo- 
geneous^ isotropic plate. The quantities (El) and (GJ ) are the 
bending and torsional stiffnesses^ respectively^ of the equivalent 
longitudinal stiffener. The procedure for determining these values 
for varying amounts of composite material in the overall stiffener 
is given in Appendix II. 

The minimum values of may be calculated from Equations (l8) and 

( 19 ) in each of the two modes of the square^, stiffened plate for 
varying values of (El)/Db ^ (GJ)/Db ^ and (A^)/bt . The critical 

stress^ ^ is then calculated from the relationship 

2 

k 7T D 

X 



The average density of the plate-stiffener system is a nondimens ional 
measure of the relative weight of the combination with the introduc- 
tion of composite materials into the system. This factor is calculated 
by averaging the densities of the system as follows 

p = P (% comp) + p, (fo base) (2l) 

*^ave *^comp base ^ ^ ' 

where comp) is the percent by volume of composite in the entire 

plate stiffener system^ and {^o base) is the percent by volume of the 
base metal used (aluminum^ steely etc.). 

The buckling efficiency of the stiffened plate combination ^^cr^^ave^ 
is plotted versus (A )/bt for varying amounts (by volume) of composite 
material in the stiffener in Figures 2 through 4 . Aluminum is the base 
material^ with either boron epoxy or graphite epoxy used as the advanced 
composite material. 

Peterson"^ determined that the optimum longitudinal stiffener for a flat 
plate has the relative proportions (l) stiffener width equal to twice 
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plate thickness^ and (2) stiffener height equal to eight times the 
plate thickness. These ratios are utilized as starting values for 
the present investigation. 
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RESULTS AND DISCUSSION 



An investigation of the axial compression buckling of a simply 
supported flat plate with a central longitudinal stiffener has been 
undertaken to estimate the effects of composite reinforcement of the 
stiffener with stiffener torsional stiffness effects included. The 
minimum or critical values of k^ ^ using Equations (l8) and (19)^ ^^e 
calculated for the two modes of a stiffened, square, homogeneous, iso- 
tropic plate for varying values of (El)/Db and plotted in Figure 5* 
These values are in excellent agreement with those of Seide and Stein^ 
for the m = 1 , n = 1 case, but the values of this investi- 
gation are slightly higher than those of Reference 5 t^ie m = 2 , 

n = 2 case. This discrepancy is due to the inclusion of the torsional 
stiffness of the stiffener in the present analysis. Without torsional 
stiffness effects, m = 2 , n = 2 case remains con- 

stant at a value of l6 despite the size of the stiffener. Timoshenko 

g 

and Gere, in the case of a plate simply supported on three edges with 
one of the unloaded edges clamped under compressive loading show the 
maximum value of 22.6 . This would be analogous to the 

present problem when the torsional stiffness is infinitely large. The 
two-term deflection function used in the present investigation will 
result in values of higher than that arrived at by Reference 8, 

but only when the stiffener size is excessively large and out of the 
area of practical interest. Equations (l8) and (l9) give values of 
which are sufficiently accurate, as the primary aim of the 
present report is to investigate the trends of plate efficiency rather 
than to calculate actual design data concerning buckling loads. 

The primary quantity for demonstrating plate buckling efficiency in this 
investigation is a /n , where cr is the critical stress of the 

plate determined from Equation (20), and is the average density 

determined from Equation (21 ). Figure 2 is a graph of cr /p for a 

square aluminum plate with a centrally located longitudinal stiffener 
augmented by varying amounts of graphite-epoxy composite. The buckling 
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efficiency has been calculated for varying stiffener size. The sharp 
knee in each of the curves results from the transition of the plate 
buckling pattern from a m = 1 ^ n = 1 configuration to a m = 2 ^ 
n = 2 mode. To the left of this knee^ or in the 1^1 mode^ sig- 
nificant increases in buckling efficiency result from the replacement 
of aluminum with graphite epoxy in the stiffener. This is due to the 
increased flexural stiffness of the stiffener^ since the graphite epoxy 
has a greater value of Young’s modulus than does aluminum. However^ to 
the right side of the knee^ in the 2^2 buckling configuration^ the 
plate buckling efficiency is actually reduced with the addition of 
the graphite epoxy. Once the 2^2 condition is reached^ the stif- 
fener lies at a node of the plate^ and the flexural stiffness of the 
stiffener has no effect on the magnitude of the critical stress. Only 
the torsional stiffness of the stiffener is a factor in this range and^ 
since the shear modulus of aluminum is approximately three times that 
of graphite epoxy^ the overall buckling efficiency decreases with the 
addition of the composite material. As an example^ a plate with an 

all-aluminum stiffener and a value of (A )/bt of 0.205 attains a 

s 

buckling efficiency parameter value of 4l.O X 10 in. (Point A). If 
the aluminum stiffener is then replaced by an all composite stiffener^ 
this parameter value is decreased to 38.8 X 10^ in. (Point B). 

It may also be noted from Figure 2 that in the 1^1 configuration^ as 
the percentage of composite in the stiffener increases^ the percent 
increase in the buckling efficiency decreases. For a plate with a 
value of (A^)/bt of O.O 8 ’^cr'^^ave 15*7 X 10^ in. for 

an all aluminum stiffener. The replacement of 10 percent of the stif- 
fener cross-section with graphite epoxy raises the value of a /p 
to 19*5 X 10 in.^ a percentage increase of almost 25 percent. Raising 
the percentage of composite to 50 percent gives a buckling efficiency 
parameter value of 28.3 X 10^ in. , or an increase of 80 percent 

over the all aluminum stiffener. An all composite stiffener has a 
/ x' 4 

u /p value of 31 •6 X 10 , or a little greater than a 100 per- 

cr' ave 

cent increase over the metal stiffener. In another case, where the 
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size of the stiffener is increased to give a value of (A^)/bt of 
0,l4 ^ a 10-percent-composite stiffener has approximately the same 
value of ^Qi-Z^ave all-composite stiffener^ although the 

latter structure will be lighter in weight and much higher in cost. 

As mentioned previously^ past studies having to do with the buckling 
properties of all-metal stiffened plates have generally neglected the 
torsional rigidity of the stiffener. Figure 5 is a graph of the same 
plate-stiffener combination as that of Figure 2; however^ the shear 
modulus of graphite epoxy is assumed to be identical to that of alu- 
minum (G = 4.0 x 10^ Ib/in. ) . In this case^ the replacement of a 
percentage of the aluminum in the stiffener with graphite epoxy always 
results in an increase in ^^rZ^ave regard to the percentage 

of composite^ size of the stiffener^ or the buckling pattern of the 
plate. Thus^ it is evident that the torsional stiffness of a rein- 
forced metal stiffener can play a significant role in the stability 
problem and should not be overlooked when high-extens ional modulus 
composites are introduced to improve structural efficiency. 

Another popular composite material in developmental use is boron epoxy. 
This composite^ although having a higher value of Young’s modulus^ is 
appreciably more dense and significantly more expensive than graphite 
epoxy. Figure 4 is a plot of the buckling efficiency of a boron-epoxy- 
augmented stiffened plate. Comparison of Figure 4 with Figure 2 reveals 
that although the use of boron epoxy results in higher values of ^cr^^ave 
than graphite epoxy in the 1,1 buckling configuration, it results in 
greater loss in efficiency when in the 2,2 mode. Higher values of 
a /p in the 1,1 mode are the result of the increased Young’s 

modulus for boron epoxy, whereas the reduced efficiency in the 2,2 
configuration is attributed to the increased weight of boron epoxy and 
the fact that the shear moduli of the two materials are essentially the 
same . 
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CONCLUDING REMARKS 



This investigation develops a buckling criterion for axially compressed^ 
simply supported^ flat^ rectangular, homogeneous, isotropic plates with 
composite-reinforced, centrally located, longitudinal stiffeners. This 
criterion is obtained from a more general development based on ortho- 
tropic plate theory. Comparison of the results of the present study 
with past investigations shows the derived buckling coefficient to be 
generally more realistic than those which neglect the torsional stiff- 
ness of the stiffener when the plate cross-section aspect ratio between 
stiffeners is at least sixty or less. 

It is also shown that the addition of composites to improve the effi- 
ciency of a stiffened flat plate requires more judicious consideration 
than simply replacing conventional material with high-performance fi- 
brous composites. Effects which may have proved to be negligible in 
the analysis of conventional structures can become significant in the 
conventional-composite combination due to the relatively large spread 
between E and G in advanced composites as compared with that in 
conventional metals. The replacement of a conventional stiffener with 
one consisting of composites will usually lead to lighter structural 
weight, but the undesirable side effects of lower buckling efficiency and 
increased cost may also be introduced. Careful study must be made as to 
the amount of composite to be used, the percentage of increase in buckling 
efficiency actually achieved, and the additional costs involved. It has 
been shown that in certain cases, a rather minimal use of composites will 
result in maximum gain in weight reduction and increase in critical stress 
with relatively small increases in cost. 

When working with advanced composites, the designer of high-performance 
vehicles can no longer neglect the torsional rigidity of a stiffener as 
was normally done with low subsonic-range all-metal systems. Such ne- 
glect can introduce errors which may tend to show that any use of com- 
posites will achieve additional gains in structural efficiency, when in 
fact this is not so. The lower shear modulus of composites requires that 
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the effects of torsional rigidity should always be taken into account^ as 
unconservative errors in critical stress may result in either premature 
buckling or reduced buckling efficiency. 

The designer must also be cognizant of all the advantages and disadvan- 
tages of the various composite materials now becoming available. As 
shown herein for graphite and boron epoxy^ each of these materials has 
certain advantages^ and the most suitable material must be chosen for 
the specific application at hand. Composite materials present the pos- 
sibility of increased strength and stiffness along with decreases in 
structural weighty but careful study must be made with regard to proper 
use and material selectivity to achieve optimized structural efficiency 
and fabr icability at minimum cost. 
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APPENDIX I 



DEVELOPMENT OF BUCKLING CRITERION 



The total potential energy, Equation (12) (with N^ = N^^ = 0 ), upon 
substitution of the deflection function, Equation (15 yields 
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Integration of this expression over the area of the plate yields 
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(25) 



Now. the first variation of Equation (25) with respect to A must 

mn 

reduce to zero since any first order change in the total potential 
must vanish for every admissible variation in the lateral displacement 
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function. Hence^ 
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Equation (24) reduces to Equation (17) when a factor of 
is removed. For the special case of an isotropic square 
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Substitution of Equation (25) into Equation (l?) yields 



A 

mn 




n a 




2(ei) m 

+ sin — 

Db 2 




P7T 

A sin — 



2(GJ) / b \ 2 2 ^ 

+ I — I n CC cos — 

Db \ri7T / 2 




bt 




ITTT 

sin — 
2 




P7T 

A sin — 
mp o 



0 



(26) 



A second approximation solution of Equation (26) for the one half-wave 
pattern in each direction yields two equations in two unknowns: namely 
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Now, for the A *s to possess other than trivial values, the deter- 
^ mn ^ ^ 

minant formed by the coefficients of the A ’s must vanish. Hence, 
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The solution of this determinant yields Equation (l8). The second- 
approximation solution of Equation (26) for the two half-wave pattern 
in each direction yields 
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It follows from the same argument as before that for a solution; 
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The solution of Equation (30) yields the result given in Equation (19)- 
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APPENDIX II 



DETERMINATION OF EQUIVALENT VALUES OF (EI) AND (GJ) 

IN COMPOSITE- REINFORCED STIFFENER 

9 

Timoshenko and Young describe a common method of transforming a com- 
posite beam into an equivalent beam of one material for the determi- 
nation of bending stiffness. This is done by multiplying the width 
of the composite section of the stiffener by the ratio of the Young’s 
moduli; that is 



E 

_ comp 
base 

where w^ is the original width of the stiffener^ and w^ is the 
modified width of the composite section of the equivalent stiffener. 

As an example, suppose ~ ^^base ' Then, as shown in Figure 6, 

the equivalent stiffener would consist of one material, the basic 
metal, but it would possess a different shape. The moment of inertia 
of the equivalent stiffener is then found and multiplied by y 

to give the equivalent value of (El) for the stiffener. 

Timoshenko and Goodier^^ give a formula for the twisting moment or 
torque of a rectangular beam in torsion as 
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However, from Reference 9^ 
torsion is known to be 
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the torsional stiffness of a beam in free 
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Substitution of Equation (52) into Equation (53) yields 
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This value is calculated for each section of the stiffener and multi- 
plied by its respective shear modulus. Addition of these components 
gives the equivalent (GJ) of the stiffener. 
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Figure 1. Geometry of Homogeneous, Isotropic Plate 
With Longitudinal Stiffener. 
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Figure 2. Buckling Efficiency Curve for Aluminum Plate With 
Longitudinal Stiffener Augmented With Varying Per- 
centage (by Volume) of Graphite Epoxy. 
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Figure 5« Buckling Efficiency Curve for Aluminum Plate With 
Longitudinal Stiffener Augmented With Varying Per- 
centage (by Volume) of Graphite Epoxy • 
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Figure 4. Buckling Efficiency Curve for Aluminum Plate With 
Longitudinal Stiffener Augmented With Varying Per- 
centage (by Volume) of Boron Epoxy. 
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Figure 5* Buckling Coefficient Curve for Aluminum Plate With 
Aluminum Longitudinal Stiffener. 
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Figure 6. Geometry of the Transformation of a Composite Beam Into 
an Equivalent Beam of One Material. 
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